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TITLE OF THE INVENTION 

CHANNEL ESTIMATION IN A CELLULAR COMMUNICATION SYSTEM 
BACKGROUND OF THE INVENTION 
Technical Field 

[0001] The present invention relates to a method and circuitry for channel 
estimation in a cellular communication system. 

[0002] The invention is particularly concerned with a cellular 
communication system in which data is transmitted as a plurality of data 
symbols over a sequence of time slots. As is known, in a CDMA system, 
data is encoded for transmission by modulating data symbols to be 
transmitted by a unique spreading code for each of a plurality of 
communication channels. Within each cell of a cellular communication 
system, spreading codes allow for a plurality of different mobile stations to 
communicate with a base station on selectively coded channels. 

[0003] When a signal is transmitted between a base station and a mobile 
station (either on the uplink or the downlink), the signal receiving unit needs 
to establish from the signal which it has received some information about the 
communication path along which the signal has travelled. This is referred to 
herein as "channel estimation" and is carried out in a channel estimation unit 
which generates a channel impulse response. Various techniques are 
known for channel estimation. The channel impulse response is required in 
order to properly decode and demodulate incoming data. 

[0004] In an article in the IEEE Journal a Selected Areas in 
Communications entitled "A Soft-Output Bidirectional Decision Feedback 
Equilization Technique for TDMA" (II (1993) September No 7) a bidirectional 
adaptive equilization technioue that takes advantage of the periodic 
transmission of training bits in the TDMA frame structure to estimate the 
location of deep fade within a time slot is proposed. 




US-A-5.732.1 14 discloses a method of detecting reference signals for a digital data 



receiver utilising a correlation process between a conjugate of a reference symbol 
and an incoming signal. 

[0005] In an earlier proposed CDMA system explained in the Technical 
Report of lEICE (1996 - 02) entitled "Variable Rate Data Transmission on 
Single Code Channels in DS/CDMA" by Okomura et al, data is transmitted at 
variable transmission rates. Channel estimation is carried out from pilot 
symbols spaced periodically in the data stream. 

[0006] One difficulty associated with this proposal is that only the energy 
from the pilot symbols is available for channel estimation, and the signal to 
noise ratio of pilot symbols tends to be poor. Thus, filtering is required to get 
more reliable estimates. In particular, in the case of mobile stations moving 
at high speeds, it may be difficult to apply a suitable filtering technique due to 
short channel coherence time. 

[0007] It is therefore advantageous to use data symbols in addition to the 
pilot symbols for channel estimation. The present invention seeks to provide 
a technique for improving channel estimation using estimated data symbols. 

BRIEF SUMMARY OF THE INVENTION 

[0008] According to one aspect of the invention there is provided a 
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method of channel estimation in a cellular communication system in which symbols 
are transmitted between a mobile station and a base station, the method comprising: 
receiving via a communication channel an incoming signal that incorporates symbols 
spread by a spreading code in a CDMA system, the incoming signal including data 
symbols and pilot symbols, a unique spreading code having been used respectively 
for pilot symbols and data symbols; 

providing a first estimate for the channel impulse response of that 
communication channel and utilising that first estimate to process the incoming 
signal to generate an estimated data symbol; 

using said estimated data symbol to generate a soft output feedback decision; 

and 

using the soft output feedback decision to make a further estimate of the 
channel impulse response. 

Th e i nvent i on i s part i cu l arly but not e xc l us i ve l y app l icab le to a CDMA system in 
wh i ch th e i ncoming s i gnal incorporat e s symbo l s spr e ad by a spr e ading cod e . 

[0010] According to another aspect of the invention there is provided a method of 
channel estimation in a cellular communication system in which symbols are 
transmitted between a mobile station and a base station, the method comprising: 

receiving via a communication channel an incoming signal which incorporates 
symbols spread by a spreading code; 

providing a first estimate for the channel impulse response of that 
communication channel and utilising that first estimate to process the incoming 
signal to generate an estimated data symbol; 
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making a further estimate of the channel impulse response by performing a 
correlation process between a conjugate version of the spreading code and the 
incoming signal, and multiplying the results of that process by a feedback decision 
based on said estimated data symbol as a weighting value for the further estimate of 
the channel impulse response. 

[0011] The incoming signal can include data symbols and pilot symbols, with a 
unique spreading code having been used respectively for pilot symbols and data 
symbols. 

[0012] According to the further aspect of the invention, the feedback decision may 
be a soft output or a hard output. 

[0013] A further aspect of the invention provides circuitry for estimating a channel 
impulse response in a cellular communication system in which symbols are 
transmitted between a mobile station and a base station, the circuitry comprising: 

a receiver for receiving via a communication channel an in coming signal 
which incorporates symbols spread by a spreading code; 

a channel impulse response generator for providing a first estimate for the 
channel impulse response of that communication channel; 

a data symbol generator for generating an estimated data symbol using the 
first estimate of the channel impulse response; 

circuitry for providing a further estimate of the channel impulse response by 

performing a correlation process between a conjugate version of the spreading code 

and the incoming signal, and multiplying the results of that process by a feedback 
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decision based on said estimated data symbol as a weighting value for tine further 
estimate of the channel impulse response. 

[0014] A still further aspect of the present invention provides circuitry for estimating 
a channel in a cellular communication system in which symbols are transmitted 
between a mobile station and a base station, the circuitry comprising: 

a receiver for receiving an incoming signal via a communication channel, the 
incoming signal incorporating svmbols spread bv a spreading code in a CDMA 
system, the incoming signal including data symbols and pilot svmbols, a unigue 
spreading code having been used respectively for pilot symbols and data svmbols: 

a channel impulse response generator for providing a first estimate for the 
channel impulse response of that communication channel; 

a data symbols generator for generating an estimated data symbol using that 
first estimate; 

a soft output feedback decision generator for using said estimated data 
symbol to generate a soft output feedback decision; 

wherein the soft output feedback decision is used to make a further estimate 
of the channel impulse response. 

[0015] The invention also contemplates a base station and/or a mobile station 
incorporating circuitry according to the first aspects of the invention. 

[0016] For a better understanding of the present invention and to show how the 
same may be carried into effect, reference will now be made by way of example to 
the accompanying drawings in which: 
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BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWINGS 



[0017] Figure 1 is a block diagram illustrating schematically symbol estimation in a 
basic CDMA system; 

[0018] Figure 2 is a block diagram illustrating a cellular communication system; 
[0019] Figure 3 is a diagram illustrating a transmission time slot structure; 
[0020] Figure 4 is a block diagram illustrating channel estimation with feedback; 
[0021] Figure 5 is a block diagram of circuitry in a channel impulse response unit; 
[0022] Figure 6 is a vector space diagram; and 

[0023] Figures 7 and 8 are simulations of acquisition probability against signal 
energy per bit to noise ratio. 

DETAILED DESCRIPTION OF THE INVENTION 

[0024] Figure 1 is a block diagram of a basic receiving system in a CDMA cellular 
communication network. Figure 2 is a block diagram illustrating a context in which 
the present invention may be used. That is, a CDMA mobile communication system 
allows a plurality of mobile stations MS1, MS2, MS3 to communication with a base 
station BTS in a common cell via RF signals on respective channels CH1, CH2, 
CH3. These channels are distinguished from one another by the use of spreading 
codes in a manner which is known perse. 

[0025] Reverting now to Figure 1, signals incoming at an antenna 30 are received 
by an RF unit 28 and supplied to an analogue-to-digital (A/D) converter 32. The A/D 
converter takes samples of the incoming signals at intervals i to generate digital 
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sampled values x. It will readily be understood that a signal may arrive at the 
receiving circuitry having experienced multi-paths with differing propagation delays. 
Each digital sampled value Xj may include therefore components of a number of 
different symbols depending on the multipath effect. 

[0026] The digital signal is supplied from the A/D converter 32 to a channel 
impulse response unit 37 which provides an estimated channel impulse response H 
for use in processing the incoming signal. A reference code generator 22 supplies a 
reference code for despreading the incoming signal both to the channel impulse 
response unit 37 and to a rake receiver 36. The reference code is denoted C and is 
the unique spreading code for the communication channel on which the signal is 
received. Thus, it matches the code which was used to spread the signal prior to 
transmission at the base station or mobile station respectively. Although Figure 1 
illustrates circuitry located at each mobile station, similar circuitry may also be 
utilised at a base station. 

[0027] The rake receiver 36 additionally receives the input digital signal x and the 
estimated channel impulse response H. It generates an estimated transmitted signal 
u which is supplied to a decoder 31. The decoder estimates symbol values s from 
the estimated signal u which has been derived from the actual incoming signal x, in a 
manner which is know per se. 

[0028] Operation of the rake receiver is known to a person skilled in the art but in 
any event is described in more detail in Application W096/24988 in the name of 
Nokia Mobile Phones Limited. 



[0029] The form in which data is transmitted is illustrated by way of example in 
Figure 3. The data symbols are transmitted in a sequence of time slots, each time 
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slot including pilot symbols (PS) and coded data. Data is transmitted as a set of 
analogue symbols which are produced from originating digital data by i,q 
modulation. In a CDMA system the symbols are spread by a spreading code prior to 
transmission. In a typical system, each time slot has a duration of .625ms. The pilot 
symbols are introduced at the beginning and end of each time slot in a coherent 
system. These symbols are readily recognisable and so the beginning and end of 
each time slot can be identified for synchronisation purposes. In the circuit of Figure 
1, the pilot symbols are additionally identified and used in the estimation of the 
channel impulse response by the channel impulse response 37. 

[0030] Figure 4 illustrates an improved version of the circuitry of Figure 1. In 
Figure 4, the estimated symbol s is fed back to the channel impulse response unit 37 
via a feedback decision generator 20 which generates a feedback decision D from 
the estimated symbol s. That feedback decision D is used to improve the estimate of 
channel impulse response generated by the CIR unit 37. An SIR block 29 receives 
the incoming signal samples x and the estimated signal u and calculates the signal 
to information ratio used to generate noise variance for the feedback decision 
generator 20. 

[0031] Figure 5 illustrates how the feedback decision D is used to improve the 
estimated channel impulse response. Figure 5 illustrates a CIR unit 37 which is 
capable of providing an estimated channel impulse response H using both fedback 
estimated data symbols and incoming pilot symbols. The circuit comprises first and 
second matched filters 10,12 which are each used to apply the relevant despreading 
code to the incoming digital signal Xin. The reference code generator 22 generates 
three codes. A spreading long code LC is applied to both of the match filters 10,12. 
A spreading short code Wp for pilot symbols is applied only to the first match filter 10 
and a spreading short code Wd for data is applied only to the second match filter 12. 
The output signals 14,16 from the first and second match filters 10,12 respectively 
are applied to an averaging circuit 18. The averaging circuit receives the decision 
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feedback D and applies that to the signals 14,16 from the matched filters to generate 
the estimated impulse response H. 

[0032] The averaging circuit 18 includes a memory 40 for holding previous channel 
impulse responses and a weighting/averaging unit 42 for taking an average of the 
earlier channel impulse responses and applying a weighting function according to 
the fedback decision value D. 

[0033] The correlation process which is used to establish the estimate of channel 
impulse response can be expressed as: 

Hn(t)= Z Dn-C* (i)x(i+t)=Dn Z C*(i)x(i+t) (Equation 1) 

i i 

[0034] In the above equation, i is the sampling time index along the received 
signal, t is the resulting sampling time index of the correlation process, x(i+t) is the 
received data sample for that index, C* is the conjugate of the spreading code and 
Dn is the nth decision feedback where n is the time index of the impulse response. It 
will be apparent that the decision feedback D may be in respect of pilot symbols 
and/or data symbols. 

[0035] In use, an initial estimate for the channel impulse response H would be 
established by the CIR unit 37 on the basis only of known incoming pilot symbols in 
a particular time slot in a manner which is known per se. Alternatively, an initial 
estimate can be made using Equation (1) without the weighting function D. Once a 
preliminary estimate of the channel impulse response has been established, this is 
held in the memory 40 and supplied to the rake receiver 36 for subsequent signal 
processing of data symbols in that time slot. The first data symbol to be estimated is 
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then subject to feedback to provide a feedback decision D so tliat a new estimate for 
the channel impulse response H can be determined using the procedure outlined 
above and as illustrated in Figure 5. The spreading code conjugate C* in equation 1 
is derived from merging the conjugates of Wd and LC (for data) or Wp and LC (for 
pilot symbols). Then, an average is taken of the new estimate and the preliminary 
estimate. Thus, at each estimated symbol an improved channel impulse response is 
obtained for subsequent processing. The average can be taken coherently within 
any coherent period, that is part of a slot, one whole slot or a multiple of slots. 
Alternatively, a non-coherent average can be taken over any period. 

[0036] It will be apparent from Equation 1 that the correlation process can 
effectively be implemented without decision feedback as if all the transmitted data 
symbols are ones. The transmission feedback D acts as an average weighting 
factor to determine the estimated channel impulse response H. 

[0037] Another aspect of the implementation discussed herein is to determine the 
nature of the decision feedback D generated by the feedback decision generator 20. 
Figure 6 illustrates the theory underlying one proposal which is to utilise a soft 
symbol decision feedback (SDFB) process by generating a soft feedback decision 
value D. Figure 6 illustrates a vector diagram of signal space. Vector A is a 
transmitted data signal. Vector A' is the received data. Vector A is the actual value 
we expect from the transmitted data signal since the rake receiver cannot be perfect. 
The received data symbol can be written as: 

A -a+n (Equation 2) 

where n is the noise vector. 
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[0038] It is natural to assume that both the signal and the noise are zero mean. 
We have: 

E{a} =0 
E{n} =0 

E{A'}=E{a+n} =0 (Equation 3) 

In a binary (+1/-1) environment, we have also: 
E{AA*} =1 

E{A'A'*} =E{(a+n)(a+n)*} 
=E{a a*}+ E{nn*} 

=E{a a*}+ (Equation 4) 

where is the noise variance. 

[0039] With these statistics of the signal, the soft output can be adjusted to satisfy 
the conditions for feedback if the noise variance is known beforehand. For a simple 
solution, we suppose that the signal loss is minimum and can be neglected: 

|a-A| 0 
Then equation (3) becomes: 
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E{AW*} =E{a a*}+ 
=E{AA*}+ 



=1 



(Equation 5) 



[0040] Using this analysis, the feedback decision generator can generate a soft 
value D satisfying the following requirements: 



[0041] As an alternative, a hard decision feedback may be utilised by performing a 
binary cut-off, that is to provide either +1 or -1 to the input of the channel impulse 
response unit for each estimated symbol S. 

[0042] The effectiveness of both solutions is illustrated in Figures 7 and 8. In 
figures 7 and 8, the x axis represents signal energy per bit to noise ratio (Eb/No) and 
the y axis represents the scale of acquisition probability - that is the likelihood that all 
signal propagation paths have been captured in the received signal as estimated 
using the C.I.R. Figure 7 is the simulation results with the best possible channel 
estimation, and Figure 8 is the simulation results in which the channel information of 
a slot is obtained only from the pilot symbols in the same slot. These simulations 
have been generated based on the MATLAB program included as Annexe 1. The 
simulation results compare the acquisition probabilities of the impulse response of 
different approaches of the decision feedback compared with no decision feedback. 
In both cases, it can be seen that the best results are given by a soft decision 
feedback for both pilot and data symbols. This is the top line in both simulations. 



E{D}-=0 



(Equation 6) 



E{DD*}= l+a^ 



(Equation 7) 



The next line is for soft feedback decisions relating only to data symbols, that is 
without pilot symbol feedback being taken into account. 

[0043] The next line is for hard data feedback decisions in respect of data symbols 
only. 

[0044] The lowest line is where there is no decision feedback. 
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ANNEXE 1 



MATLAB program 

% Soft Decision Feedback 

% ARIB specification 

% QZD/NTC - FEB. 98 - Oulu 

clear, rand ( ' seed' , 0) , randn ( ' seed' , 0) 

Dat_Weight = 0.5; SNR = -30:3:0; Fingers = 10; 

h=[0.7+j*0.1,0,0,0,-0.5-j*0.3,0,0,0,0,0.3+j*0.5,0,0, ... 

0.1-j*0.13,0,0,0,0,0,0,0.3-j*0.1,0,0,0,0.5+j*0.5, ... 

0,0,0,-0.5-j*0.3,0,0,0,0,0.3+j*0.5,0,0,0.1-j*0.13, . . . 

0,0,0,0,0,0,0.3-j*0,l] ; 
h=h/sqrt (h*h' ) ; 

AcqError = zeros (size (SNR, 2 ) , 4); 

for kkk = l:size (SNR,2) 

snr = SNR (kkk) % SNR in dB 

for nnn=l:1000 
DAT = (rand(l,40)>0.5)*2-l; 
% CODE 

CL_I=(rand(l,2560)>0.5) *2-l; CL_Q= (rand ( 1 , 2560 ) >0 . 5 ) *2-l ; 
C = CL_I + j*CL_Q; 

w=l;for k=l:6;w=[w w;w -w] ; end; wO=w ( 1 , : ) ; wl=w (2 , : ) ; clear w; 
W0=[];W1=[]; for k=l:40, W0= [WO , wO ] ; Wl= [Wl , wl ] ; end 
CO = wO .*C; CI = Wl .*C; 
% Spreading 

D=[]; for k=l:40, D- [D, DAT (k) *wl] ;end % Data 
P=ones (size (D) ) ; % PILOT 

xl = (D + j*P) .* C; 
% Channel 

x2 = conv(xl,h) /2; 
% Noise 

nA=l/10^ (snr/20) ; f or k=l:size (x2,2) 

n (k) =nA*randn*exp ( j *rand*2*pi) ;end 
x = x2 + n; 

%%%%%%%%%%%%%%%%%%% MF %%%%%%%%%%%%%%%%%%% 
IRPilot = zeros (39, 64); IRData = zeros (39, 64); 
for k=l:39, l=(k-l)*64; 
for m = 1 : 64 

IRPilot(k,m) = IRPilot (k,m) + x (l+m-l+ (1 : 64) ) *C0 (1+ (1 : 64) ) ' 
IRData(k,m) = IRData(k,m) + x (l+m-l+ (1 : 64) ) *C1 (1+ (1 : 64) ) ' 
end 

end 

IRPilot = -j*IRPilot/64; 
IRPilotHard=sum (IRPilot (1:28, :) ) /28; 
IRData = IRData/64; 

%%%%%%%%%%%%%%%%%%% H' %%%%%%%%%%%%%%%%%%%%%%%%%% 
a = IRPilotHard ,* conj (IRPilotHard) ; 
H = zeros (size (h) ) ; 
for k = 1: Fingers, findf = 0; 
while ( findf < 1) 

[n,m] =max (a) ; a (m) =0; 
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if (m<=size (h,2) ) , H (m) =IRPilotHard (m) ; findf=2; end 

end 

end 

%%%%%%%%%%%%%%%%%%% decision %%%%%%%%%%%%%%%%%%% 

hl=find(abs (h)>0.01) ; 

h2=hl; % Idea 

%h2=find(abs (H) >0. 01) ; % blind 

PilotSoft = zeros (1,28); DataSoft = zeros (1,40); 

for n=l: Fingers 

a = x(h2 (n) :size (C0,2)+h2 (n) -1) conj (CO) * conj (H (h2 (n) ) ) ; 
for k=l:28 

PilotSoft (k) = PilotSoft (k) + imag (sum(a ( (k-1) *64+l : k*64) ) ) ; 
end 

a = x(h2(n) :size(Cl,2)+h2(n)-l) conj (CI) * con j (H (h2 (n) ) ) ; 
for k=l:40 

DataSoft (k) = DataSoft (k) + real (sum(a ( (k-1) *64+l : k*64) ) ) ; 
end 

end 

PilotSoft = PilotSoft/64; 
DataSoft = DataSoft/64; 
DataHard = sign (DataSoft) ; 

%%%%%%%%%%%%%%%%%%% Averaging %%%%%%%%%%%%%%%%%%% 
a = 0; b = 0; c = 0; 

IRHard=zeros (1, 64) ; IRSoft=zeros (1, 64) ; IRSof tp=zeros ( 1 , 64 ) ; 
IRSofter=zeros (1, 64) ; IRSof terp=zeros (1, 64) ; 
for k=l:39 

IRHard=IRHard+IRData (k, : ) *DataHard(k) ; a=a+abs ( DataHard (k) ) ; 
IRSoft=IRSoft+IRData (k, : ) *DataSoft (k) ; b=b+abs (DataSoft (k) ) ; 

end 

for k=l:28 

IRSoftp=IRSoftp+IRPilot (k, : ) *PilotSof t ( k) ; c=c+abs 
(PilotSoft (k) ) ; 

end 

IRSoftp = IRSoft*Dat_Weight/b + IRSoftp * ( l-Dat_Weight ) /c; 
IRHard = IRHard*Dat_Weight/a + IRPilotHard* (l-Dat_Weight) ; 
IRSoft = IRSoft*Dat_Weight/b + IRPilotHard* ( l-Dat_Weight ) ; 

probability 
a = IRPilotHard .* con j ( IRPilotHard) ; 
for k=l:size(hl,2) , [n,m] =max (a) ;h2 (k) =m;a (m) =0;end 
c=0; for k=l:size (hl,2) , for n=l : size (hi, 2) 

if (h2(k)==hl(n) ) , c=c+l; h2(k)=0; end 
end, end 
AErrNoDFB = size (hi, 2 ) -c; 
a = IRHard conj (IRHard) ; 

for k=l:size(hl,2) , [n,m] =max (a) ;h2 (k) =m;a (m) =0; end 
c=0; for k=l:size (hl,2) , for n=l : size (hi, 2) 

if (h2(k)==hl(n) ) , c=c+l; h2(k)=0; end 
end, end 
AErrHardDFB = size (hi , 2 ) -c; 
a = IRSoft .* conj (IRSoft) ; 

for k=l:size(hl,2) , [n,m] =inax (a) ;h2 (k) =m; a (m) =0; end 
c=0; for k=l:size (hl,2) , for n=l : size (hi, 2) 
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if (hi (k)==hl (n) ) , c=c+l; h2(k)=0; end 
end; end 
AErrSoftDFB = size (hi , 2 ) -c; 
a = IRSoftp .* conj (IRSoftp) ; 

for k=l : size (hi ^ 2 ) , [n,m] =max (a) ; h2 ( k) =m; a (m) =0; end 
c=0; for k=l : size (hi, 2) , for n=l : size (hi, 2) 

if (h2 (k)==hl (n) ) , c=c+l; h2(k)=0; end 
end, end 
AErrSoftPDFB = size (hi , 2 ) -c; 
AcqError (kkk, : ) =AcqError (kkk, :)+.,. 
[AErrNoDFB,AErrHardDFB, AErrSoftDFB, AErrSoftPDFB] 
end % nnn 
end % kkk 
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